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Abstract 

Tow-dimensional diffusers are used for a variety of applications, such as turbo machines. 
An analytical procedure for computing the static and total pressure distribution subsonic diffuser 
with axi-symmetric or two dimensional cross-sections used .Finite element procedures for 
incompressible viscous fluid motion, turbulent, unsteady flow are used in this paper. The unsteady 
flow problem is solved by the Newton-Raphson method. Straight-walled, two-dimensional 
diffusers for several aspect ratio and angles were investigated theoretically for the purpose of 
studying the regime of no appreciable stall. Introducing the (k-s) model to the work also used the 
effect of turbulence. 

1. Introduction: - 

The deceleration of flowing fluid, or diffusion as it is widely known, is a process of 

primary importance whereby some of the kinetic energy of the flowing fluid is converted into static 
pressure rise. Hence, fluid motion take places against an adverse pressure gradient, and unless great 
care is taken in the design of diffusing system. From the point of view of the fundamental of fluid 
motion, the study of diffusing flows is indeed very interesting. Therefore, it has attracted the 
attention of many investigators, and a great deal has been published on performance of various 
diffuser configurations, e.g., conical, annular, rectangular, radial. 

There is a wealth of literature on flow through the diffuser. Many researchers presented 
design and analysis of the diffuser, both experimental and analytical solution. Riyadh [1], presented 
the mathematical model for the two-dimensional diffuser in centrifugal compressor. The flows 
described in this research steady state, and turbulent flow. K. Eisele, A. Ongoren,, P. Holbein and 
M.V. Casey [2], described the dedicated test rig for planar diffuser flows of different opening angles 
is available and has been modified to take the flow conditions such as inlet blockage and inlet 
turbulence level of a typical turbo -machinery diffuser into account. The flow in this diffuser was 
analyzed with LDA at different opening angles in the range of attached and stalled flow regimes. 
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The measurements show that a three-dimensional separation of the flow sets-in a corner at moderate 
opening angles and grows with increasing diffuser angles to a big separation zone downstream. Kai 
U. Zieglerl, Heinz E. Gallus and Reinhard Niehuis [3], presented the interaction between impeller 
and diffuser is considered to have strong influence on the flow in highly loaded centrifugal 
compressors. However, the knowledge about this influence is still not satisfying. This two-part 
paper presents an experimental investigation of the effect of impeller-diffuser interaction on the 
unsteady and the time averaged flow configuration in impeller and diffuser and the performance of 
these components. The flat wedge vaned diffuser of the investigated stage allows an independent 
adjustment of diffuser vane angle and radial gap between impeller exit and diffuser vane inlet. 
Attention is mainly directed to the radial gap, as it determines the intensity of the impeller-diffuser 
interaction 

The diffusing passage two dimensional diffusers considered here is show in Fig. (1). 
Geometry is fixed by aspect ratio, [b/Wl], and any two of the for dimensionless geometric 
parameters, (20, L/Wl). The computational procedure taken solves the time average Naviar-Stokes 
equation for unsteady state, incompressible and turbulent flow in two-dimensional situations. 
Analysis of the flow inside the diffuser was based on finite element technique, by using Galerkin 
method. 




Fig. (1) Two-dimensional diffuser for centrifugal compressor. 



2. Design and parameter of diffuser: - 

The simple two dimensional diffuser may be described by three geometric parameters, the 
area ratio (AR), the divergence angle of one way (0), and the ratio of the wall length to entry width 
(LAV). The parameters are illustrated in Fig. (1). But only two of the three will be independent since 
they are related by Ref. [4] as: 

AR = l-2 — Sin 3 .... (1) 

W 

It has been common practice to plot of diffuser characteristics against (20) and (L/W). An 
alternate method would be to use (AR-1) and (L/W). Several combinations have been used in this 
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paper, since each has its own advantages. 

The geometric of simple vaned island type diffuser system consists of a number of channels, 
as shown in Fig. (2). The channel can be divided into three regions, the flow region (1) is unsteady 
and very unsettled because of jet and wakes emanating from rotating impeller passage and these 
tend to mix as the flow advances towards the diffuser throat. Region (2) clearly conforms to parallel 
wall rectangular diffuser. In region (3), as in region (1), the flow is bounded by the diffuser blade 
only on one side. Hence, the mechanism of diffusion would not be the same as in the region (2), as 
described by Ref. [5]. 

Because the distance from the wall is very important, therefore takes path (1), and it's at 
(0.02m) distance from the wall. In this paper, it is intended to design region (2) only, in designing 
the diffuser, several models under line (A- A) of Fig. (3) Were taken, as illustrated by [5, 6 and 7]. 




Fig. (2) Wedge type vane Island diffusers Fig. (3)Flow Regimes for two-dimension 

diffuser. 

3. Variation Equations 

The basic conservation equations are the continuity and motion which for unsteady state, 
take the following time averaged form in two-dimensional case [8 and 9]. 

■Continuity 

*ij=0 -..(2) 

■motion 

p[^ + u ij u jj- (7 ijj= 0 -(3) 

Where: 

tr.. = (—PS- + lUffd:;) and d i; = —(u, ,+«..) 

y N y r e Jj y ' 2 

The well known Navier-Stokes equation can be obtained by introducing equations ( cr ) and 
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( d tj ) into equation (3) and using equation (2). 

( du } 

P {^ + U iJ U jj + P i- M eff U iJj = ° - (4) 
■ Two equation turbulence modeling 

In this paper the turbulent effects are included via a k-s model, which is based on the 
transport of two additional turbulent quantities, the turbulent kinetic energy (k) and its dissipation 
rate (s) for high Reynolds number, with two dimensional unsteady flow. These two additional 
balance equations may be written in the following form, [8 and 10] as: - 

(duf_ 



DK d 

P = — 

Dt dy 



eff dk 



+ M eff{dy \ -** - (5) 



Ds d 
P — = — 
Dt dy 



k ■ j 



eff ds 
a E dy 



(du) r pe* 



-C.-M- ^ -C 0 ^— ...(6) 
U effydy) 2 k 



Where 



(ldf=fll+flT and |u T =Cnpk 2 /8 

At high Reynolds numbers, 8 may be assumed proportional to (pk 3/2 /f) Here, Ok, Ci, a £j C2 
and Cu are the turbulence model constants, which are respectively equal to 1.0, 1.44, 1.3, 1.92 and 
0.09 respectively. 

• Pressure Recovery Coefficient: - 

This is the most widely used parameter and it gives static pressure rise as a function of the 
inlet velocity head, i.e.; 

CP = ^T ...(7) 

2 t 

4. Finite Element Analysis 

To apply the finite element method, variation equation corresponding to equation (2) and (3 
are required. They can be obtained by employing the conventional procedure called Galerkin 

method. Let u* be the weighting function. Multiplying both side of equation (3) by u* , 

integrating over the whole volume V and using Greens theorem , the following is obtained, [9, 11 
and 12]. 

J p — l -u- dV = \u- Uj -u -dV + \u- <t» -dV =Q . . . (8) 

v dt v ,J J v J,J 

Introducing equations ( cr. ) and ( d tj ) and rearranging them, the final form of the variation 
equation is derived. 

\p — l -u-dV = \ pu-u- -u -dV - \(Pu. .)dV + u^\u. . \u. . +u . . \dV = 0 •••(9) 
As to the equation of continuity, after multiplying both sides of equation (2), by weighting 
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function and integrating over the whole volume V, the variation equation is expressed in the 
following: 

\{P*u- -)dV = 0 ... (10) 

v l 'J 

Assume that the flow field to be analysis is divided into small regions called finite 
elements, let the interpolating equation for velocity and pressure in each finite element be expressed 
by the following form, [12]: 

P = Va ...(12) 
Where O a and ^ denote interpolation function or shape function for velocity and pressure 
respectively. ( ) means the velocity at ath nodes of each finite element in the ith direction, and 

( P x ) is the pressure at the Xth node. Using the finite element mesh, the flow inside diffuser has been 
studies. 

The boundary conditions are given as follows Fig. (2). 

u=0 and v=0 for all nodes at (AB and CD) . ...(13) 

u=100 and v=0 for all nodes at AC. ... (14) 

P=0 for all nodes at BD. ... (15) 

For weighting function , and P- , the relation that is similar to equations (11) and (12) are 

used as follows. 

u* =Q> a u *. ...(16) 

l CCl 

* * 

P ...(17) 

Introducing equations (11), (12), (16) and (17) into equations (9) and (10), the following 
finite element equations can be obtained for motion and continuity as: 

M cd^^ +K apf^ +H cdA + S cd(3j u pj= Q •■ < 18) 

H M u a =0 ...(19) 



where 



And superposed dot denotes partial differentiation with respect to time (t). The conventional 
finite element superposition procedure gives the final equation system for the whole flow field in 
the following form. 

M n v*R + K n v 0 v +H 0 q 0 +S fl v o -0 ...(20) 
ap r apy (5 y aA A ap (5 v ' 

H a ,v a =0 ...(21) 

where: - 
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v a and qx denote the velocity and pressure of the whole nodes in the flow field. 
By the same method solve the Kinetic energy equation and dissipation equation for the two 
dimensional flow equations (5 and 6) 

The type of the element in the mesh is triangular with (144) element and three nodes for each 
element Fig. (4). 



Triangular 
Element 




X 



Fig. (4) Natural coordinate triangular element. 



5. Computer Program: - 

A general purpose computer program Finite element mechanics is to producing the results to 
be discussed. The program is written and developed in Quick Basic language to solve the two 
dimensional unsteady flow Ref. [11]. The perturbation method is employed to solve the unsteady 
flow problem by equations (20 and 21). Assume that the term Q a can be expanded into Taylor 
series in small perturbation parameter (£) as follows. 

Expanding velocity vp and pressure q^ as: 

v j3= v p +£vl j3 +£2v } - (23) 

q j3 =q ( ] 3 +£q l j3 +£ 2 qj l ... (24) 

Introducing equations (22, 23 and 24) into equations (20 and 21) and equating the 
coefficients of the same order terms in (£) lead to the following simultaneous equation system : 



6 



For the nth order 



ap V apy J3 y a/C 1 ^ ap ft 



(25) 
(26) 



Replace the differential with respect to time t by the difference with respect to the short time 
increment At as: 



v fi- It 



(27) 



Where v^(0) denotes the initial value of the velocity increment for short time increment At. 

Introducing (27) in (26) get the matrix and solve this matrix by Newton-Raphson method to get on 
velocity component and pressure at all nodes. 

6. Results and Discussion 

Several models were taken for different angles (9), area ratios (AR) and the ratio 
length/width (L/W, ), the models taken are in table below: 



Model 
No. 


AR 


L1AV1 


L2AV1 


L3AV1 


L4AV1 


91 


92 


93 


94 


1 


1.4 


3 


4 


5 


6 


7.63° 


5.73° 


4.58° 


3.82° 


2 


1.6 


3 


4 


5 


6 


11.42° 


8.58° 


6.87° 


5.73° 


3 


1.8 


3 


4 


5 


6 


15.2° 


11.42° 


9.15° 


7.63° 


4 


2 


3 


4 


5 


6 


18.9° 


14.25° 


11.42° 


9.53° 


5 


3 


3 


4 


5 


6 


36.86° 


28.07° 


22.62° 


18.92° | 


6 


4 


3 


4 


5 


6 


53.13° 


41.11° 


33.4° 


28.07° | 



Fig. (5) Presented the velocity component at x-direction for path (1). The velocity decreases 
with increase (X/L) because of increased in divergent angle for the same angle and the velocity 
increases with decreased angle for the same (X/L). 

Fig. (6) Shows the velocity component at y-direction along the diffuser for different angles 
for path(l) and area ratio (1.4). The velocity decreases along the diffuser for the same angle due to 
increase in divergent duct. And the velocity increases with decreased the diffuser angle for the 
distance between (0.1-0.5) and then no effect with increase or decrease angle. 



Fig. (7) Presents the pressure recovery coefficient, distribution along the diffuser for path 
(1). The pressure recovery coefficient increases with the increased distance for the same angle due 
to the increase in the static pressure because of the decreasing velocity (diffusion process). Pressure 
recovery coefficient decreases with the increased angle due to the increased absolute velocity and 
the decreased static pressure and it began from negative value near the wall due to high velocity at 



7 



inlet. 

Fig. (8) Shows the velocity component at x-direction with area ratio, for path (1) and 
different (L/Wl) The velocity component decreases with increase area ratio for the same (X/L) due 
to increase in the divergent angle . And the velocity component increase with increase (X/L) for the 
same area ratio specially when the area ratio above (2). 

Fig. (9) Presented the velocity component at y-direction with variable area ratio for path (1) 
and different (X/L). The velocity component increased slowly with increased area ratio due to 
increased in divergent angle for the same (X/L). And velocity component decreased with increased 
(X/L). 

Fig. (10) Shows the pressure recovery factor with different area ratio. The pressure recovery 
factor increased with increased area ratio at the same (X/L) due to increased in divergent angle. And 
it increased with increased (X/L) at the same area ratio due to increased in angle. 

Fig. (11) Shows area ratio with divergent angle (20) for Path (1) and different (X/L). The 
area ratio increased with increased divergent angle for the same (X/L). And area ratio increased 
with increased (X/L). 

Fig. (12) Present velocity component at x-direction with area ratio for three paths (1,2 and 
3) and for (X/L =3). The velocity component decreased with increased area ratio for the same path. 
And the velocity component decreased with increased distance from the wall at the same area ratio. 

Fig. (13) Fig. (12) Present velocity component at y-direction with area ratio for three paths 
(1,2 and 3) and for (X/L =3). The velocity component increased with increased area ratio for the 
paths (land 2) and then decreased due to increased in divergent angle. And the velocity component 
approaches to constant and equal zero for path (3) due to increased distance from the wall and 
increased velocity component at x-direction. 

A compression between the numerical results from the present mathematical model with 
other workers analyses was made. Fig. (14), Fig. (15) and Fig. (16) Shows a compression between 
the results obtained with that published by Riyadh [1]. 

Fig. (14) Shows the velocity component ratio at y-direction for diffuser for AR=2, L/W1=6 
and path (1). Noted that it have different between the present model (unsteady) and Riyadh [1] 

d 

(steady) due to effect of an unsteady term ( — ). 

dt 

Fig. (15) Present the effect of pressure recovery coefficient across the diffuser for AR=2, 
L/W1=6 and path (1). It has been found the effect of pressure recovery in steady state is smaller 
than unsteady state. 

Fig. (16) Shows the velocity ratio for AR=2, L/W1=6 and path (1). The velocity ratio 
decreased with increased distance due to increased in divergent angle but when unsteady state, the 
velocity decreased less than steady state due to effect the term of unsteady. 
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10 20 30 40 50 

20 

Fig.(ll) Effect of area ratio for variable 
divergent angle of diffuser for path(l). 




AR 

Fig.(12) Effect of velocity component at x-direction 
for variable area ratio and varaiable paths. 





AR=2....Path(l) 

- Present (unsteady) 

- Riyadh (steady) , 



1.0 1.5 2.0 2.5 3.0 3.5 4.0 

AR 

Fig.(13) Effect of velocity component at y-direction 
for variable area ratio and varaiable paths. 



X/L 

Fig.(14) The velocity component at y-direction 
for two-dimensional diffuser for path (1). 




AR=2....Path (1) x 

- Present (Unsteady) 

- Riyadh (Steady) i 



l — i — i — I — i — i — i — r 



0.0 



I -1 

0.8 



AR=2....Path(l) 
- Present (Unsteady) 
Riyadh (Steady) , 




0.2 0.4 0.6 

X/L 

Fig.(15) Effect of pressure recovery coefficient 
accros diffuser for path(l). 



1 I 1 1 

0.2 



n -1 

0.4 



X/L 



Fig.(16). The velocity component at x-direction 
for two-dimensional diffuser for path (1). 
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7. Conclusion 



This paper has dealt with a finite element method to analyze unsteady flows of viscous fluid in 
two dimension of diffuser. Discretization procedure is based on Galerkin approach.. The following 
conclusion can be made: - 

1 . From studying the results of analysis of flow throw the diffuser for several models, the geometry 
data for optimum design found are: 

AR=2, LAV, =6 and 0 = 9.53° 

2. This model gives a pressure recovery coefficient (CP) equal to (0.75) and low velocity 
component (V) with respect with other models and lows ratio OF (LAV1) with other models. This 
agrees well with that published by reference [1,2, and 3]. 

4- Pressure recovery coefficient increase with increase the angle of divergence of diffuser wall and 
increase in static pressure with decreased in average velocity. 



Nomenclature 



Ci, C2, C K :- Constant of turbulence. 
K: -kinetic energy of turbulence. 
I :-Length scale of turbulence, (m) 
Pi>Inlet static pressure. (N/m 2 ) 
P2:-Exit static pressure. (N/m 2 ) 
Cp:-Pressure recovery coefficient. 
L: -Length of diffuser. 
v:-Kinematics viscosity. (m 2 /s) 
u:- velocity component, (m/s) 
t:-Time (sec.) 

9:-Divergent angle. (degree) 
Q:-Area of element, (m 2 ) 
a,X,p,y:-The global node number. 



a:-Stress tensor. 

Wl and W2:-Inlet exit diffuser width between diverging wall, (m) 

b:-distance between parallel walls of a diffuser. (m) 

a4SS:- Diffuser vane suction side angle (Fig. 1) 
r 2 , r 4 :- Radius of impeller. 



(j,e:-Laminar viscosity. (kg/m.s) 
|u, t :- Turbulent viscosity. (kg/m.s) 
s:-Dissipation rate of turbulence. 
\i e ff: -Effective viscosity. (kg/m.s) 
<3k:- Turbulence Prandtle number for k-equation 
a s : Turbulence Prandtle number for s-equation 
R:-Radial component. 
Vi=Absolute velocity 

Cd=coefficient equal to unity without losses. 
X:-axial distance, (m) 
p: -Density. (Kg/m 3 ) 
x:-Shear stress. (N/m 2 ) 
i, j:-Tensor notation subscripts. 
AR: area ratio. 
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